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Expressions for the Wick contractions contributing to the scalar pion form-factor were computed 
model-independently in chiral perturbation theory at next-to-leading order The results reveal 
correlations amongst the different contractions in terms of low-energy constants and allow for 
extrapolating lattice data for individual Wick contractions. The quark disconnected contribution 
to the real part of the form factor turns out to be suppressed with respect to the quark connected 
one. The corresponding contribution to the scalar radius has the same size as the connected 
contribution and can therefore not be neglected. 
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1. Introduction 

Quark disconnected Wick contractions constitute a considerable numerical problem for lattice 
QCD. An estimate of the full quark-propagator is needed and the signal-to-noise ratio typically 
leaves a lot to be desired. In current simulations in the iso-spin symmetric limit disconnected con- 
tractions contribute to a number of phenomenologically relevant quantities (g — 2, ;r7r-scattering, 
hadronic /T-decays, . . . ). The situation becomes worse once also iso-spin breaking effects are in- 
cluded in the simulations since disconnected contributions will then contribute in a larger variety 
and to many more quantities. 

Given the difficulties in computing quark disconnected contractions to a satisfactory precision 
in lattice simulations one often neglects them, thereby introducing an unknown systematic effect. 
Large efforts therefore go into devising dedicated algorithms for improving the numerical evalu- 
ation of quark-disconnected diagrams (e.g. ^ ^). Here, an analytical and model-independent 
approach for computing individual Wick contractions in chiral perturbation theory which was pre- 
sented in [||] is applied to the case of the scalar pion form factor for QCD with A'^^ = 2 and 
Nf = 2 + 1 dynamical flavours. At next-to-leading order (NLO) in the chiral expansion the con- 
tribution of the disconnected contraction to the form-factor is numerically small. Both the discon- 
nected and the connected contribution however contribute with roughly the same magnitude to the 
scalar radius. The connected and disconnected contributions turn out to be coixelated in terms of 
low-energy constants of the chiral Lagrangian. The expressions derived here can be used to guide 
chiral extrapolations of lattice results for individual Wick contractions. As an aside the argument 
allowing to compute the quark-connected part with partially twisted boundary conditions [^, ||] for 
an improved momentum resolution for the scalar form factor is presented. 

2. Technique 

The scalar form factor of the pion is defined as 

{7i'{p')\uu + dd\Kjip)) = 5'jFsait), (2.1) 

where t = {p' — p)^ is the squared momentum transfer between the initial and final pion and the 
sub-script on the r.h.s. identifies the scalar form factor for Nf = 2 flavour QCD. In lattice QCD the 
matrix element on the l.h.s. is computed in terms of the ground-state contribution to the Fourier 
transfomi of the Euclidean 3-pt. function 

(^0'{z)S{y)0'\x)) , (2.2) 

constructed of the interpolating operators 0'{x) = Yi^'jsYiix) and S{x) = \j/2{x)Y2{x)- The = 
{u,d) are SU (2) flavour vectors of u- and J-quarks and the matrices T,- = a,/2 are proportional to 
the Pauli matrices. Following the two types of Wick contractions contributing to this coiTclation 
function are 

{uysddddjiu) = [Uy^vvddyiu) + {uy^dvvdysu) . (2.3) 

Note the additional valence quark v with m„ = = that was introduced in order to construct 
the decomposition on the r.h.s.. At this stage the connected contribution consists entirely of flavour 
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off-diagonal currents and the ai^gument presented in Q allows to compute at least this contribu- 
tion using partially twisted boundary conditions, thus improving the momentum resolution. Each 
individual term on the r.h.s. represents an unphysical correlation function in the unphysical theory, 
QCD with an additional valence quark v. The expression in chiral effective theory for the ground 
state contribution to the l.h.s. has been derived many years ago in chiral perturbation theory at 
NLO [0 ^ and later at NNLO [|l|, ^ in both SU{2) and SU{3) chiral perturbation theory. The 
effective theory frame work for the computation of the corresponding contributions to the r.h.s. is 
partially quenched chiral perturbation theory (PQ;tPT) [T5|] which is a frame-work allowing 
to vary the sea-quark and valence-quark content independently. While the individual terms on the 
r.h.s. represent unphysical matrix elements, their sum represents a physical process. It is therefore 
conceivable to use different methods for the computation of the terms on the r.h.s.. The connected 
contribution is conveniently computed in a lattice simulation. Where the disconnected contribution 
is difficult to compute numerically the method advocated here can be appUed. 



3. PQ;^ PT for the scalar form factor 



According to [|14|, [ISp the chiral effective theory for Nf = 2 and Nf = 2 + I QCD with an 
additional valence quark v is developed around the graded flavour symmetry groups SU (3 1 1 ) and 
SU (4| 1), respectively. The leading order chiral Lagrangian is [|^ [14, 15, 16, 17] 



(3.1) 



where X = 2B{s + M). The mass matrix has the form M = diag {mq,mq,mq,mq) in SU{3\1) and 
M = dia.g{mq,mq,ms,mq,mq) in SU{4\1) and we define the external scalar- source as 5 = 2T"s" 
(with T" a generator of the flavour group). The relevant counter terms can be derived from the 
Lagrangian, 



= L4Str{d^,U{dt'Uy}str[xU^ + Ux^} +L5Str[{d^U{d''Uy){xU'^ + Ux' 
+L6Str\^XU^ + Ux^y + UStr^^Ux^Ux^ + ZU^XU^f ■ 



(3.2) 



The remaining calculation proceeds as in standard chiral perturbation theory. The known result for 
the full form factor is reproduced. 



Fl^it) = 2B 



+ 



(2f- ml] 



B{ml,t) 



(3.3) 



Fl,{t) = 2b{ 1 + ^2 ( - + + A3 + + ,0 + '-B{m\,t)) } , 

The expressions for A„,2 and B{m^,t) are standar^d and can be found in [^. The new results ar^e the 
individual expression for the connected (F^j^^) and disconnected contribution (F^j^^): 



Fi:2(0=2B{l + i,(-iA„,. 
=2b{0 + -1( 



,D , (f + ml) , 



K'nl,t))} : 



(3.4) 
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if3(f) = 2b{i + -!( - + iA,„2 +A? + + '-B{ml,t) + 'I^B{m^,mlt)) } , 

=2B{0 + i,( +A^+^I±^B{ml,t) + '^B{m'^,t)^'fB(ml,ml^^^^ 

Note that for both Nf = 2 and Nf = 2 + I, = + F^, as expected. The A^^^ contain combi- 
nations of the low energy constants (LECs), 

Al =4{m2(-8L:j-4L^+16L^+8L^)+?(2L^ + L0}, 

4|m2(-4L^-4L^+8L^ +8L0+? L^}, (3.5) 

4{m2(-4L^ +8L^) +? 2^ }, 

Af = 4{m2(-6L^-4L^+12L^+8L§)+m|.(-4L^ + 8L^)+f(2L^ + L0} , 
AC = 4{m|(-2L^ - 4L5+4L^ +SLl)+ml{-4Ll + 8L^)+f } , (3.6) 

A? = 4{m2(-4Li +8L^) +f 2L^ }. 



AC 



A^^ 



and 



The are the LECs of SU (3) chiral perturbation theory [|16| and the L- are related to the better- 
known SU (2) LECs /[ via matching, 

Al=tll +4n4/5, 

AC = tir^-SLD+Amlir^ +4Ll-SL'^) , (3.7) 

A^ = t{ +8L^)+4m|( -4L^+8L^). 

The remaining SU (3|1) constants and Lg ai^e less well known since the corresponding terms in 
the L(4)-La grangian can be removed in the case of the flavour group SU (2) (trace-identities). For 
Nf = 2+1 similar expressions for the individual Wick contractions can be derived for the octet and 
singlet form form factors 

{7i'\uu + dd-2ss\n'') = 

r3 8^ 

{7i'\uu + dd + ss\7i'') =d'''F^{t). 



From the expressions in eq. ( |3.4| ) we see that while the connected contribution to the form 
factor starts at leading order, the disconnected contribution starts at NLO and is therefore expected 
to be smaller in magnitude. Turning to the scalar radius (r^) = 6^^^^\t=Q, the derivative with 
respect to the momentum removes the leading term from the connected contribution and therefore 
connected and disconnected contribution start contributing at the same order in the chiral power 
counting. 

Concentrating on the Nf = 2 + 1-case we now fix the free pai^ameters in the expression for 
F^l^^^ {t) with input as summarised in table |l|. The results for the 2-flavour case, which turn out 
to be very similar both qualitatively and quantitatively, can be obtained after fixing the LJ" through 
matching the unphysical LECs to the 3-flavour theory [H]. The results are shown in the plots in 
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Table 1: Values for the 3-flavour low energy constants used in the illustration of the results. The subtraction 
scale is ;U = 0.77GeV. 



figure |l] for the full form factor and the connected and disconnected contribution, respectively. 
For the real part the magnitude of the disconnected contribution is sub-dominant in line with the 
above observations. Above the two-pion threshold the disconnected contribution also contributes 
to the imaginary part and is similar in magnitude to the connected contiibution. Since the leading 
contribution has no imaginary part this does not come as a surprise - the imaginaiy pait in both the 
connected and the disconnected contribution start at the same order, NLO. We also show the mass- 
dependence of the scalar radius. What could have been anticipated from the first of the three plots 
by observing that the slope around ? = of the connected and disconnected contribution is very 
similar, is confirmed in the third plot - the connected and disconnected contributions to the scalar 
radius are nearly of the same size. For a lattice simulation this underlines that the disconnected 
contribution cannot be neglected [|l|, §]. 



4. Conclusion 



Partially quenched chiral perturbation theory proves to be a powerful tool for understanding 
in detail the dynamics in the low-energy sector of QCD. The decomposition of correlation func- 
tions into individual Wick contractions provides guidance for lattice computations in various ways: 
First of all an estimate of the magnitude of the disconnected contribution can be made in a model- 
independent way. Secondly, NLO LECs can be extracted from only the connected contribution to 
the form factor which is numerically much better accessible in lattice simulations. The expression 
does however contain new linear combinations of LECs and one has to study whether the LEC one 
is interested in can be extracted. Thirdly, a full computation of the form factor in lattice QCD is 
still very challenging due to the numerical cost of computing the disconnected contiibution to a 
satisfactory precision. In this situation one can either compute only the connected contribution in 
lattice QCD and predict the disconnected one, or one computes the connected contiibution for a 
large set of parameters, while the disconnected one only for a reduced set of parameters to fix the 
LECs and to eventually extrapolate it to the physical point. 
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Figure 1: First row: the plots show the momentum dependence of the real (left) and imaginary (right) parts 
of the scalar form factor for Nf = 2 + I: Full form factor (solid red), connected contribution (dotted blue) 
and disconnected contribution (dashed green). The bottom plot shows the scalar radius for Nf = 2+1 as 
a function of the squared pion mass (niK is fixed to the physical value and ~ |m|. — ^m^, same color 
coding as first row). 
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